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Polymer melts without topological constraints are be- 
lieved to be described by the Rouse model [ij, i.e. the 
random forces from the bath acting on the center-of-mass 
(CM) rjq{t) of a reference chain of length N are supposed 
to be uncorrelated. For the corresponding CM mean- 
square displacement (MSD) hj^{t) this implies 
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for all times t with d being the spatial dimension, sa 
b'^W/N the self-diffusion coefficient and b the effective 
bond length pj]. The monomer mobility W may be ob- 
tained by fitting the monomer MSD 

h{t) = {{r{t) - r(0))2) = b^{Wt)^"' with a = 1/4 (2) 

(r(t) being the monomer position) for times smaller than 
the relaxation time Tn ~ N'^/W 1]. In fact, correlated 
(colored) forces have been observed both experimentally 
and numerically as reviewed in Ref. [2]. These colored 
forces are best captured in a computer simulation using 
the "velocity correlation function" (VCF) Q 
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St° for t > St, 
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which measures directly the curvature of /iN(i) with re- 
spect to time with u{t) = r^{t + St) —r^{t) being the CM 
displacement for a time window St. Using simple scaling 
arguments it has recently been suggested 3^ that an in- 
terplay of melt incompressibility and chain connectivity 
implies a negative algebraic decay of the VCF 
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for i < Tn « iVVW^ (4) 



with c being an empirical amplitude, p the monomer den- 
sity and io — (2 + d)a the exponent characterizing the 
colored forces. This scaling has been demonstrated nu- 
merically by means of Monte Carlo (MC) simulations for 
three-dimensional (3D) melts where w = 5/4 Q. The 
aim of this Note is to show that Eq. (U) holds also in 
effectively two-dimensional (2D) polymer films. 

As in Ref. we use a version of the bond-fluctuation 
model (BFM) with finite excluded volume penalty e and 
topology violating local moves to the 26 neighboring lat- 
tice sites. Length scales are given in units of the lattice 
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FIG. 1: Various MSDs for polymer melts confined to ultrathin 
films of width H — 4 computed using the BFM algorithm with 
topology violating local MC moves The open symbols 
represent the monomer MSD h{t) which compare well to the 
Rouse prediction (dashed line). Systematic deviations from 
the Rouse model are revealed for the CM MSD h-M{t) (filled 
symbols) as emphasized by the effective power-law exponent 
/? = 0.9 (bold line). The stars correspond to Eq. ((5| with 
c = 1 and = 20N^ for N = 8192. 



constant, time in Monte Carlo Steps (MCS). The chains 
are confined to a slit of width H — A between two par- 
allel repulsive walls following Ref. Q • All properties dis- 
cussed refer to their 2D projection parallel to the walls. 
Working at the 3D volume fraction (f> = 0.5 of occupied 
lattice sites [H, the relevant 2D number density is thus 
p = {cj)/8)H = 0.25. The systems are effectively c? = 2 di- 
mensional since the width of the slit is much smaller than 
the typical chain size. For the overlap penalty £ = 10 (in 
units of k^T) presented here the static properties are sim- 
ilar to Ref. [J] where £ = oo. Note that the monomers can 
still cross each other in the direction perpendicular to the 
walls. The chains are thus "self-avoiding trails" [aL i.e. 
they adopt logarithmically swollen configurations [J, [1] ■ 
Since R'^/N diverges for large N (albeit very weakly), the 
operational definition of the effective 2D bond length b is 
a delicate issue and th e pr ojected effective bond length 
of the bulk, b = 3.24-\/2/3 ~ 2.65, is used as reference. 
Note that the scaling of the dynamical properties dis- 
cussed below does not depend on this specific value. 

It is important to stress that due to the use of finite 
overlap penalty and topology non-conserving MC moves, 
the dynamics remains to leading order of Rouse-type as 
can be seen from the MSDs presented in Fig. [T] The 
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FIG. 2: Scaling plot of the directly computed displacement 
correlation function CN{t) for a large range of chain lengths A'^ 
(open symbols) using the parameters p, b, W and Dn ~ 
indicated in the plot. The data collapse for t ^ Tn confirms 
that CN{t) ~ The corresponding time exponent cj = 1 

(bold line) is observed over up to five orders of magnitude in 
time. The exponent uj = 2 for overdamped colloids expected 
for t ^ Tn is indicated by the dashed line. 

monomer MSD h{t) (open symbols) compares well over 
several orders of magnitude to Eq. 1^ indicated by the 
dashed line. Although it is likely that this scaling does 
not hold in a strict sense due to the mentioned logarith- 
mic swelling of the chains, deviations — if they exist — 
are apparently of higher order and irrelevant on the log- 
arithmic scales we focus on. As can be seen for N — 64, 
the monomers diffuse again freely (thin solid line) for 
times t 3> 7n. Please note that it was not the aim of the 
present work to sample for our larger chains (TV > 2000) 
over the huge times needed to make this regime acces- 
sible. The short-time behavior, Eq. can be used to 
determine W w 0.002 for the effective local mobility and 
to predict a self-diffusion coefficient Dn ~ 0.005/iV us- 
ing the Rouse model 3 . Systematic deviations from the 
Rouse scaling are revealed, however, for the CM MSD 
h-[^{t) (filled symbols) as emphasized by the effective ex- 
ponent /3 = 0.9 (bold line). At variance to the 3D bulk 
[3| deviations are visible up to a time t* scaling as the 



chain relaxation time Tn. (It is thus only possible to con- 
firm the given value of Djq from the long-time plateau of 
h^{t)/2dt for chains up to TV = 1024.) 

A more precise characterization of the colored forces 
acting on the chains is achieved by means of the VCF 
Cn (t) computed directly as described in Q . The rescaled 
data presented in Fig. [2] confirm the predicted algebraic 
decay, Eq. ^ , with an exponent u ~ 1 and an empirical 
amplitude c « 1 (bold line). Note that the scaling im- 
plies CN(t) ~ i-e. the correlations are due to ^ 
local independent events as for 3D melts. The large time 
behavior for t ^ Tn can be seen for the shorter chains 
which approach the exponent CL' = (2-|-d)/2 = 2 (dashed 
line) expected for overdamped colloids in incompressible 
solutions Using Eq. (0) it follows from Eq. (g]) that 

h^{t)^AD^t(^ + ^^\n{t*/t)^ fort<TN (5) 

with t* being an integration constant. Please note that 
the logarithmic contribution to Eq. ([S]) dominates for t <^ 
t*. While in 3D we have t* - < Tn it turns out 
that for thin films t* « 2QN'^ sa Tn/18 yields a good fit 
for all chain lengths we have sampled as can be seen for 
N = 8192 in Fig. [1] (stars). 

In summary, we have presented here computational 
work on the displacement correlations in thin polymer 
films without topological constraints and momentum 
conservation using a well-known lattice MC algorithm 
d, 0| . As anticipated by the scaling result suggested re- 
cently for 3D melts [Sj] , it is shown here for effectively 2D 
systems that the VCF Cn (t) reveals a negative algebraic 
decay with a power-law exponent w = (2-|-d)/4 = 1. This 
implies a logarithmic correction to the CM MSD. Our 
MC approach corresponding to a strongly overdamped 
melt should be relevant to most experimental setups due 
to the strong frictional forces from the walls. Obviously, 
qualitatively different correlations are to be expected if 
hydrodynamic interactions matter — as is the case for 
molecular dynamics simulations of a bead-spring model 
with a weak Langevin thermostat Q — or if the chains 
are confined to strictly 2D layers due to the ensuing com- 
pactness and surface fractality of the chains [H, |^ . 
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